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p • Abstract 

The off-shell dynamics of the 0(3) nonlinear sigma-model is probed in terms 
^ . of spectral densities and two-point functions by means of the form factor 

■ approach. The exact form factors of the Spin field, Noether-current, EM- 

tensor and the topological charge density are computed up to 6-particles. The 
corresponding n < 6 particle spectral densities are used to compute the two- 
point functions, and are argued to deviate at most a few per mille from the 
exact answer in the entire energy range below 10'^ in units of the mass gap. 
To cover yet higher energies we propose an extrapolation scheme to arbitrary 
particle numbers based on a novel scaling hypothesis for the spectral densities. 
It yields candidate results for the exact two-point functions at all energy scales 
and allows us to exactly determine the values of two, previously unknown, 
non-perturbative constants. 
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1. Introduction 



The O(A^) nonlinear sigma (NLS) models describe the dynamics of 'Spin' fields 5* = 
{Si, . . . , Sn) taking values in the — 1 dimensional unit sphere and governed by the 
action 



where go is a dimensionless coupling constant. Classically it has two important sym- 
metries: First the invariance under the action of the internal O(A^) group rotating the 
spins, and second the invariance under spacetime conformal transformations. The QFT 
is thought to describe an 0(A^)-multiplet of stable massive particles, the mass scale being 
non-perturbative in the coupling constant [|l|, ^. Thus, the first of the above classical 
symmetries is preserved, but the second is lost. 

It may be worthwhile to briefiy recapitulate the physical picture behind this phenomenon, 
which is most succinctly done in the lattice formulation. The action (1.1) is replaced by 
a discretized version on some finite L x L square lattice. The functional integral is then 
well-defined and can be approximately evaluated by means of Monte Carlo simulations. 
The continuum limit is taken by driving the system into a critical point. For the bare 
coupling constant it is generally believed that the only critical point is go 0. (See 
however p, |^ for the possibility of a Kosterlitz-Thouless (KT) type phase transition 
for go > 0.) There exist non-trivial spin configurations whose energy goes to zero as 
L (yo for fixed go and which therefore in infinite volume are present at arbitrarily 
small go- These so-called super- instanton configurations can be thought of being the 
"enforcers" of the Mermin- Wagner theorem ^. They disorder the spins, forbidding 
a spontaneous magnetization in two dimensions. In QFT language this means that the 
O(A^) symmetry is unbroken and (in the absence of a KT-phase transition) the theory has 
a mass gap. It also means that perturbation theory (PT), known to be renormalizable 
in finite volume 0, has infrared problems as it starts from a fictitious ordered ground 
state. Nevertheless for O(A^) invariant correlation functions a result by Elitzur and David 
guarantees that (using periodic or Dirichlet boundary conditions) the coefficients of 
the bare lattice PT expansion have finite L — oo limits. The bare PT expansion can 
then be converted into a renormalized one, and the renormalized one into an expansion 
in the running coupling constant, where the running is defined through the perturbative 
beta function It must be emphasized that all this can be (and is) done regardless 
whether or not the final expansion coefficients have any relation to, and significance for, 
the unknown exact correlation functions. The hypothesis of asymptotic freedom is that 




(1.1) 
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they do have. Namely, if one were given the exact correlation functions and tried to 
perform an asymptotic expansion in the running coupling constant (the running again 
defined by the perturbative beta function) the claim is: (i) this expansion exists and (ii) 
the coefficients obtained coincide with the ones computed in PT. 

This is a mathematically precise statement which is either true or false. It is true that it is 
commonly believed to be true - though not everyone might wish to take this as a substitute 
for a proof. Indeed, the correctness of the claims (i) and (ii) has been challenged in a 
series of papers P, §], stimulating some controversy Evidently the problem cannot 
be addressed within the perturbative frame. It is also difficult to even come close to 
conclusive results on the basis of Monte Carlo simulations alone; the presently available 
data cover only the energy range below 50m, where m is the mass gap of the model. This 
is far too low to study the elusive high energy properties. The purpose of this and an 
accompanying technical paper is to study the off-shell dynamics of the 0(3) NLS model 
via the form factor approach. Form factors in this context are the matrix elements of some 
local operator between the physical vacuum and some asymptotic multi-particle states. 
In a QFT with a factorized scattering theory they can in principle be determined exactly 
by solving a set of recursive functional equations || that use the exact two-particle S- 
matrix as an input. Once the form factors are known, off-shell quantities, e.g. two-point 
functions, can be computed by inserting a resolution of the identity in terms of multi- 
particle states. For models whose S-matrix is diagonal in isospin space, the resulting low 
energy expansion has been seen to converge rather rapidly |Tl], |12|. Models with a 
non-diagonal S-matrix are technically much more demanding. However they also provide 
a much better testing ground for 4-dimensional QFT scenarios and deserve more dedicated 
attempts to gain insight into their exact off-shell dynamics. 

In the 0(3) NLS model we computed the exact form factors of the Spin field, Noether- 
current, EM-tensor and the topological charge density up to 6-particles. Their two-point 
functions are evaluated by means of a Kallen-Lehmann spectral representation [|T3]. The 
n < 6 particle spectral densities are argued to provide results for the two-point functions 
that differ only about 3 to 10 per mille from the exact answer in the entire energy range 
below lO^m — lO^m, depending on the quantity considered |14|. In the low energy range 
the agreement with the Monte Carlo data of Patrascioiu and Seller [15| is excellent. At 
higher energies we compare with renormalization group improved 2-loop PT. Within the 
range considered 2-loop PT yields an (within 1%) accurate description of the system only 
for energies above 50m — 100m, provided one uses the known exact value of the Lambda 



parameter [|l^] to fix its absolute normalization. Based on the low energy Monte Carlo 
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data alone, however, one would be tempted to maximize the apparent domain of validity 
of PT by tuning the Lambda parameter such as to match the relevant part of the Monte 
Carlo data. Doing this in the NLS model the Lambda parameter comes out wrong by 
about 10%. Generally speaking this emphasizes the importance to have an independent 
estimate for the onset of the (2-loop) perturbative regime. In the case at hand this is 
provided by the form factor results. 

A major challenge remains the computation of the extreme UV properties of the model 
independent of PT. For the important case of a two-point function this amounts to sum- 
ming up all multi-particle contributions to their spectral resolution. The short-distance 
asymptotics of the two-point functions is related to the /i ^ oo asymptotics of the cor- 
responding spectral density p(/i). We approach the problem by taking advantage of a 
remarkable self-similarity property of the n-particle spectral densities: For large n and 
A > 1 they appear to behave like 

P^^^H/^) - ^ WH^/'"" P'"-^ (mWm)V^^"") , (L2) 

where m is the mass gap and 7, a are certain 'critical' exponents. Given p^^\ji) for some 
initial particle number n, (1.2) can be used - without actually doing the computation 
(which for large is practically impossible) - to anticipate the structure of all > n 
particle spectral densities. We promoted (1.2) to a working hypothesis and explored its 
consequences. The results obtained and the status of the hypothesis may be surveyed as 
follows: 

• It is a statement about the spectral densities of the exact theory at all energy scales. 

• Its consequences for the UV behavior are consistent with PT. 

• It produces new exact non-perturbative results like the normalization of the spin 
two-point function at short distances. 

• It allows one to compute numerically the two-point functions at all energy/length 
scales. 

The paper is organized in the following way. In section 2 we describe general features 
of the spectral representation of two-point functions and their asymptotic expansions. 
Applied to the 0(3) NLS model we prepare various results on the four operators under 
consideration. The next section contains a summary of our results on the 0(3) form 



3 



factors, form factor squares and their properties, as well as an exact expression for the 
asymptotics of the ri-particle spectral densities. The results on the two-point functions 
based on the n < 6 particle form factors are described in section 4. The extrapolation to 
arbitrary particle numbers by means of the above scaling hypothesis is implemented in 
section 5, to be followed by brief conclusions. 



2. Spectral representation of two-point functions 
2.1 Spectral representation 

The form factors characterize an (integrable as well as non-integrable) QFT in a similar 
way as the ra-point functions do. Assuming the existence of a resolution of the identity 
in terms of asymptotic multi-particle states, the ra-point functions can in principle be 
recovered from the form factors. In the important case of the two-point functions this 
amounts to the well-known spectral representation. For the Minkowski two-point function 
(Wightman function) of some local operator O one obtains in a first step 

W^'ix -y) = J2l-ff[^ g-^(xO-,0)p(")(e)-.(.i_,i)p(")(e) |_^(n)(^)|2 ^ ^2.1) 

where = (0|O(0)|^„, ...,61) are the form factors of O and P^TKO) = 'Z^P^l{0^), 

with Pq{9) = mchO, pi{6) = msh.6 are the eigenvalues of energy and momentum on 
an n-particle state.Q The local operators are classified by various quantum numbers, 
in particular by their Lorentz spin s and their mass dimension A. It turns out to be 
convenient to parametrize their form factors in terms of 'scalarized form factors', which 
are functions of the rapidity differences only and carry quantum numbers A = s = 0. For 
an operator with quantum numbers (A, s) we shall use the parametrization 

J^(")(^) = >C(p("n^)) /^"H^) > (2-2) 
^Our kinematical conventions are: {x^,x^) are coordinates on 2-dimensional Minkowski space TR^'^ 
with hilmear form x-y ~ x^'rii^iyy", 77 = diag(l, — 1). Lightcone coordinates are — {x° ±x^)/V^ = Xz^z] 
the norm is ||a;|| = y/±x ■ x = V ±2x+a;~, ±a;^ > 0. The antisymmetric tensor is eoi = — cio = 1- 
The normahzation of the 1-particle states is {9i\02) — An5{9i — 02 ), which corresponds to the standard 
normahzation in d + 1 dimensions, speciahzed to c? = 1. For simphcity we assume that aU particles are 
of the same mass m and suppress internal indices in section 2.1. 
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where in the cases of interest £ is a polynomial of degree A > |s| in the n-particle momenta 
P^"^(^) of integer spin s, i.e. 

£(p(")(^)) = e-""/:(p(")(0 + u)) . (2.3) 

The Wightman function (|2.1|), considered as a distribution, can then be obtained by 



iy^(x -y) = C{td,) C{idy) W{x - y) , (2.4) 



differentiation 

W''{x-y) = C{td,)C{iOyj 

where W{x-y) is defined as the r.h.s. of (^) with |J^(")(^)p replaced by |/(")(^)p. We 
shall refer to W{x) as the "scalarized Wightman function of O" . Let us emphasize that in 
general W{x — y) can not be interpreted as the two-point function of some local (scalar) 
field; it is only a useful auxiliary function from which the physical 2-point function can 
be obtained through differentiation. 

For many purposes it is useful to rewrite (|2.1j ) in terms of a Kallen-Lehmann spectral 
representation [|T^. Changing integration variables according to 

'm(e^^ + . . . + e^")\ 



Ui = 9i — 9i^i , l<2<r;, — 1, a; = ln 



u] = m 



M^^){u) 

n + 2 ^ ch(Mj + . . . + Uj^i) 



(2.5) 



one obtains for the scalarized Wightman function 

roo 

W{x-y) = -i dfx p{fx) D{x-y;fx) , 
Jo 

P(/i) = E P^"H/^) ' p«(/i) = 5(/i-m)|/«r 



n>l 



p^-Hp) = f '^''\^;;^tr' i/^"^^)^ ^(^ - m^-\u)) , n>2. (2.6) 

Notice that no problem of convergence arises for the spectral density. First, each n- 
particle contribution exists because, for fixed /i, the integrand has support only in a 
compact domain V{p) C 1R"~^ in which the form factors are bounded functions so that 
the integration is well-defined. Viewed as a function of yU one observes that p^'^\jj) has 
support only for p > mn. Summing up the n-particle contributions to p{fi) therefore only 
a finite number of terms (those with n < [p/m], [x] being the integer part oix) contribute. 
Under some mild assumptions on the growth of p(/i) (for example it is sufficient to require 
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that p(/i) is polynomially bounded in /i) the existence of the spectral density guarantees 
that of the 2-point function, as defined through ( p.6| ). The integration kernel is given by 



D{x; m) = ] 9ix') [sign(xo) Mm\\x\\) - iYo{m\\x\\)] + ^e{-x^)Ko{m\\x\\) , (2.7) 

and coincides with the 2-point function of a free scalar field of mass m. We use the conven- 
tions of [1^ for the Bessel functions. The differentiation (2]^) is a bit cumbersome in the 



general case; usually however one will be interested in the behavior at spacelike distances 
in which case only higher order modified Bessel function K„(m||x||) arise. Alternatively 
one can use the Fourier representation of D{x; m) and write 

W'^ix -y) = dfipift) I ^ 9{po)S{p' - /i^) dp) C{-p) e-^ (--^) . (2.8) 

On general grounds the Fourier transform of W'~^{x) must have support only inside the 
forward lightcone = {p E IR^'^ | > 0, po > 0}. From (|2.8|) one finds indeed 

wo,, /^£(P)£(-P) peV* 

w^{p) = < WpW 

[ p^V+ 

Thus, up to kinematical factors the spectral density can also be viewed as the Fourier 
transform of the two-point Wightman function. 

For comparison with perturbation theory one needs the time-ordered two-point function 
and its Fourier transform. Its spectral representation is easily read off from ( |2.6|) , ( |2.8| ) 



G^(x - y) = I ^ e-'^<--y^ C{p) C{-p) n-p' - le) , (2.9) 



where 



/(.)= ^d^^^. (2.10) 
Jo z + jJ'^ 

The definition of I{z) has been chosen such that it has a cut along the negative real axis 
and one can recover the spectral density from the discontinuity along this cut 

p(/i) = — disc /(— /i^) . 
vr 

Conversely, disc I{z) determines I{z) up to a polynomial ambiguity. 

The spectral representation of the Euclidean two-point function (Schwinger function) is 
obtained similarly. The Schwinger function can be defined by S'^{xi,X2) ■= W'-^{—ix2, Xi) 
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for X2 > and then by analytic continuation to X2 < 0. By construction it also coincides 
with the analytic continuation of G'^(x). In a momentum space integral this formally 
amounts to the replacement {po,pi) = {ipi ,pf)- The Euchdean counterparts of ( p.8|) and 



(|2.9| ) thus are 

S^ix -y) = / ^ e-^^-(^-^) C^p) C^i-p) I{p') 

= C''{id,)C'^{idy)S{x-y) . (2.11) 

We shall also use the notation {0{x)0{y)) for 5*^(2; — y). The integration in (|2.11| ) is 
now over the Euclidean momenta (pi,P2) = (pf^pf) and C^{pi,p2) ■= C{ip2,Pi)- The 
"scalarized Schwinger function of O" entering the second line is 

Six) = i- rf/ip(/x) Ko {^i^xl + xl) . (2.12) 

Notice that the spectral density p(yu) and hence the function I{z) is the same in the 
Minkowski and in the Euclidean situation. Given either p(/i) or I{z) for a specific oper- 
ator all two-point functions can be computed - scalarized and physical, Minkowski and 
Euclidean ones. 

As an example consider the case of the energy momentum (EM) tensor. It is of additional 
interest, because its spectral density is closely related to the Zamolodchikov C-function. 
In 1 -|- 1 dimensions, the symmetric and conserved EM tensor can always be parametrized 
in terms of a non-local scalar field r, the "EM-potential" 

T^y{x) = C^u{idx) t{x) , C^y{p) := -p^,py + r]^^p^ . (2.13) 

The defining relation ( p. 2]) for the scalarized EM form factors thus reads 

{o\T^M\dn,.--,o,) = [-p^r\o)P!)''\o)+r]^,p^-\e)p^-^'ie)] f-\e^...e,) , (2.14) 

and the scalarized form factors can be interpreted as that of the scalar field r. Using ( p.6|) 
the spectral representation of the Minkowski two-point function is 

POO 

{Q\Tf,^{x)Tc,p{y)\Q) = -iC^y{idx)Ca(3{idy) j d^i p{fi) D{x - p) , (2-15) 

and similarly for the time-ordered two-point function. Their analytic continuations to 
xo = —1x2 coincide and yield the Schwinger function according to the rules ( |2.11|) , (|2.12|) 
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with = -PaPis + ^app^. Explicitly 

{T,,{x)T^M) = I (^^""'^'"'^ iP^P-^ - V P'\ [PcPp - Sa(^ P'] I{P') 

= /:^Md^)C^^{zdy)S{x-y). (2.16) 

The lightcone components correspond to the 5*0 (1, l)-irreducible pieces. In particular 
and the trace t = 2T^ have polynomials = — and 2£_|__ = p"^ respectively. 



Combining ( p.7|) and ( p. 61) one finds for the behavior at small spacelike distances 



(0|T++(x)T++(0)|0) = ^ (^^) (^ + ---' ^'<°' ^2.17) 

where the dots stand for less singular terms whose form is mo del- dependent. The coeffi- 
cient c is given by 

0=1271 dfipifi) (2.18) 
Jo 

and coincides with the central charge of the Virasoro algebra in the conformal field theory 
describing the UV fixed point of the model. This latter fact is part of the statement of 
Zamolodchikov's C-theorem The proof of the C-theorem is particularly transparent 

from the viewpoint of the spectral representation in the Euclidean case [|19[. In particular 
using (|2.16|) , ( p.l2|) , the properties of the modified Bessel functions and the definition 



( |2.18| ) it is easy to arrive at the "C-theorem sum rule" |^| for theories with a mass gap 



c = Stt y d^xx^ (t(x)t(O)) 
where {t{x)t{0)) is the Euclidean correlator of the trace. 



2.2 Asymptotic expansions 

Suppose that in a QFT the validity of asymptotically free perturbation theory has been 
established in the sense outlined in the introduction. Then PT allows one to compute 
the asymptotic expansion of the Fourier transform of the two-point functions for large 
Euclidean momenta. This can be translated into the language of the spectral density 
using the defining equation (|2.10|) . Motivated by the asymptotically free case, we define 

^In particular c is finite; for other operators the integral over the spectral density will in general not 
converge. 
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^(0 by 



z \ m 

and R{0 by 



I{z) = -il\n^] (2.19) 



p(/i) = -i? An^V (2.20) 
/i \ mJ 

Using the relation ( p.lO| ) we can compute the asymptotic expansion of R{C) in terms of 



and its derivatives: 

oo 

i?(0^E^2s+i^('^+^^(0, (2.21) 
where Hi = 1, = —^({2) and all the numerical coefficients H2S+1 are polynomials in 

00 -1 

ak) = Y.-, k>2. (2.22) 

n=l 



In (|2.21|) the symbol ^ signifies that both sides of the equation have the same asymptotic 
expansion in 1/^, i.e. their difference is smaller than any power of 1/^ for large ^. In 
addition to the asymptotic relation ( p.21| ), there is a relation between the integral of -R(0 



and the leading term in the asymptotic expansion of (provided, of course, both are 
finite): 

/■oo /"OO 

l{oo)= d^R{^)= dfipifi). (2.23) 
JO Jo 

In perturbation theory the asymptotic expansion is usually given with the help of some 
running coupling constant, in terms of which the expansion is a power series and which 
goes to zero as the momentum goes to infinity. There are many functions having these 
properties. A simple and convenient choice is e(^) defined by 

^ +«:lne(0 = e, (2.24) 



e(0 



where the parameter k is related to the first two beta function coefficients of the model. 
We shall refer to e(^) as the universal running coupling function because only the universal 
part of the beta function enters its definition. Any other running coupling function can 
then be expressed as a power series in e(^). In particular this would hold for the lattice- 
measurable running coupling (for various possible definitions of the latter, see |pT|). 

If the expansion of the function l starts off at the p^^ power of e 

6(0 ^ioeP + ii e^+' + i2 e^+' + • • • (2.25) 
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using ( |2.2lD we can write 



R{0 ^ ro eP+^ + ri 6^+^ + 6^+^ + ■ • • , (2.26) 

where 

n = -pKLo- {p + 1)li, 

r2 = ^ioC(2)p(p+l)(p + 2)-pK^o-(p+l)«:^i-(p + 2)i2. (2.27) 



2.3 0(3) Nonlinear Sigma-Model 

The 0(3) non-hnear cr-model is asymptotically free in perturbation theory and admits 
an exact on-shell solution via the the bootstrap method. The asymptotic single particle 
states are assumed to consist of an 0(3) triplet of massive particles with mass m: \a,6), 
a = 1,2,3. In this paper we will study the properties of the four most interesting local 
operators in the 0(3) model: The spin field the Noether current J^, the EM tensor 
Tfj_^ and the topological change density q. In the following we shall discuss these four 
operators consecutively. 

Spin field: We normalize our spin field operator $'*(x) by 

(0|$"(0)|6,^) = (2.28) 

is proportional to the Lagrangian field renormalized in (say) the MS scheme 

<^%x) = C S^{x) , (2.29) 

where ( is some unknown finite renormalization constant. The n-particle form factor of 
is defined by 

(0|$'^(0) |a„. Or,;...; a,, 6,) = ^ (^„, ...,9,), {n odd) . (2.30) 

Here we introduced the normalization factor 2/y^ for later convenience (c.f. section 
3.2). The subscript '1' on the r.h.s. indicates that under an 0(3) rotation these functions 
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are intertwiners S*^" 3 onto the irreducible representation 3 of isospin 1. Only the 
odd particle form factors are non-vanishing because ^'^ is odd under the internal parity 
reflection of the 0(3) variable. In addition they depend only on the rapidity differences. 
As a consequence one can write 



6^'Ft\u):= /iV-a.(^n,...,^i)7t„.(^n,...,^i). (2.31) 

ai...a 



The spectral density is given by 



p^P''^(/.) = -^pr+^)(/i), (2.32a) 



^ fc=0 



pi"^(/^)= r '^";4;,w!r"' sii. - m^-\u)) . (2.32b) 



/o (471)*^- 1 
From p^P™(/i) or its Stieltjes transform 



' z + jjfi 



pp-^^{z) = I dfi " , ^7 (2.33) 



all two-point functions of the Spin field can be computed. 
For example the Schwinger function is 

" I (1^ e-^P-(^-^) rP^'^(p2) . (2.34) 
In the time-ordered Minkowski two-point function —iPP^'^[—p'^ — ie) would enter. The 



normalization (|2.28|) is thus equivalent to the usual (re-)normalization condition on the 



Fourier transform of the time ordered Minkowski two-point function 

- « /^P'°(-p2 - ie) = + . . . , p2 ^ . (2.35) 

p2 _ jj^Z _|_ 

It also implies that for large (spacelike) distances all two-point functions decay like in the 
free theory. In particular for the Schwinger function this means 

1 



S'^'^x) = -y== e-^' [1 + 0(l/r)] , r = ^xl + xl. (2.36) 
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We will later determine the spin-spin spectral density and two-point function by the form 
factor bootstrap method. We will then be able to compare their asymptotic behavior 
with the results obtained in perturbation theory. We shall denote the universal running 
couphng function for the 0(3) model by a{p), which is defined as in the previous subsection 
to be the solution of 

^ ■ Ina(p) = In — . (2.37) 



a{p) m 

(The K parameter of (|2.24|) is equal to unity for the 0(3) model.) Perturbation theory at 



2-loop order predicts for the large p asymptotics 



p2 jspm(p2) ^ + (2 + eo) + (2 + eo)«(p) + O {a\p)) I . (2.38) 



Here the overall constant Ai is related to the unknown finite renormalization in ( p.29| ) 



while the parameter gives the connection between the perturbative mass parameter 



Aj;jg and the exact mass gap m. In the 0(3) model the latter is known exactly |jT6| 



Tfl 

^0 = In-— = ln8- 1 ^ 1.07944. (2.39) 



Finally, using the results ( |2.25D , ( |2.26D and ( |2.27| ) we obtain the perturbative prediction 



for the asymptotic expansion of the spin-spin spectral density 

p^P'^(/.) = - {l + a(/x) + O ia\i,))} . (2.40) 

In section 5 we shall find the value of the undetermined non-perturbative constant Ai to 
be 

A. = A . (2,41, 

Note that the perturbative series for both the Fourier transform and the spectral density 
contain integer powers of the running coupling only. This fact, which will play an impor- 
tant role in what follows, is a special feature of the 0(3) model, since for the otherwise 
rather similar O(A^) models with > 3 one has 

p2ppi-(p2) _ [a(p)]-iv^ p''''"(Ai) ~ - [a(Ai)]^ • (2.42) 

A* 

Noether current: In terms of the Lagrangian variable the Noether current is 

ja ^ ]_^abcgbQ^gc _ (2.43) 

9o 
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In the QFT we fix the normahzation such that its time components satisfy the equal time 
commutation relations 



JSiO, x) , Jo^(0, y)\ = le-"'^ ^,{0, y) 5(x - y) . (2.44) 

In 1 + 1 dimensions a conserved current can always be parametrized in terms of a (non- 
local) scalar field, the "current potential" j"" via 

J^,=e^.ad"r- (2.45) 

The scalarized form factors of can be identified with that of the field j". Their defining 
relation is 

(0| j;(0)|a„, On,...; a,, 6,) = C,{P^''\e)) /i%„...,,(ft„, • • • , ^i) , (n even) (2.46) 

where C^{p) = —ie^^p'^- Note that no confusion can arise here from using the same symbol 
fian --ai in (|2.3CI|) since there it is defined for an odd number of particles only. We also 
note that ( p.44|) fixes the residue of the two-particle form factor 

<2paa2ai 

fla.a. (02, = + . . . , e^^e. + ZTl. (2.47) 

(72 — (71 — tTT 



If we now extend the definitions ( p.31| ), ( p.32[ b) for even values of n as well, the spectral 



density in the case of the Noether current is given by 



oo 



k=l 

Finally, the Fourier transform for the current is defined as in (|2.33|) with the operator 
superscript '^p™ replaced by Using p™"^"^ and I™"^"^ one can represent all the two-point 
functions of the Noether current. For example the Schwinger function is 

72 

Mb /" " P -ip-(x-y) r _r)2A 1 n^l 



(r.i^) Jliy)) = ^ / ^"^"'^'"'^ [pmP. - r i'''"{pn . (2.49) 



Perturbation theory predicts for the Fourier transform [23 



13 



and, using ( |2.25| ), ( |2.26| ) and ( ^.271 ), for the spectral density 

P'^'^'if^) = ^ {l + + O («'(/i))} . (2.51) 

Topological charge density: In terms of the Lagrangian variable the topological charge 
density reads 

q = ^e"''' e^"" S" d.S'' d.S" . (2.52) 

OTT 

We shall use a different symbol for its Euclidean version 



1 



k = ^e"'' eg) 5" d^S' dpS' , (2.53) 

(where = —e^^ = 1) to emphasize that it differs by an extra factor —i from the 
Euclidean continuation of g: k = —iq^^\ This extra factor of —i, which is a consequence 
of the linearity of q in the time derivative, is responsible for the important fact that the 
Euclidean two-point function of is a strictly negative function for all non-zero Euclidean 
distances. 

A result of the form factor approach is that the topological charge density operator can 
be parametrized in terms of a dimensionless non-local scalar field $ as follows (c.f. section 
3) 

^=(□ + ^2)$. (2.54) 

The scalarized form factors of q can be identified with that of $. Their defining relation 
is 



(0|g(0)|a„, . . . ; ai, ^i) = ^oCiP^^^^O)) Ua^-aM, • • • , ^i) , {n odd) (2.55) 

where C{p) = —p^ + m? . The subscript '0' indicates that under an 0(3) rotation these 
functions are interwtiners 3®" 1 onto the singlet (isospin zero) irreducible represen- 
tation. The normalization of the topological charge operator between physical states of 
fixed particle number is not known a priori. In section 5 we shall determine the non- 
perturbative normalization constant to be 

Ao = ^ . (2.56) 
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Introducing the squares 



Ft\u)= E l/0,a„...a,(^n,...,^l)|' (2.57) 
ai...a 



we can write the spectral density for this operator as 

oo 

p*°P(/i)=AoEpf^'H/^), (2.58a) 



A:=0 



Finally if we define similarly as in (|2.33| ), the Euclidean two-point function of the 



Euclidean topological charge density is given as 

{k{x) k{y)) = - / e-^-(--^) (p2 + ^2)2 jtop(p2^ _ (2.59) 
Perturbation theory predicts 

/ip'^-^l/i) = Y^p(/i) + 0(«3(/.))} , (2.60) 

where 4°'') being the coefficient of a contact term in the two-point function, cannot be 
calculated in PT. 

Energy-momentum tensor: Since we already discussed the EM tensor in section 2.1, here 
we can be brief. From the Lagrangian one obtains 



T,, = - [d^S'^d^S'' - - r,^, . (2.61) 

In the QFT we fix the normalization such that the Hamiltonian ^'^^dxT^Q^x) has eigen- 
value ^/p^ + on an asymptotic single particle state of momentum p. The scalarized 
form factors can be identified with that of the EM-potential r as introduced in (|2.13| ). To 
fix the notation we repeat their defining relation 

(0|T^,(0)|a„, On,...- ai, ^i) = C^,{P^''\e)) fo,a„-aAOn, . . . , ^i) , {n even) (2.62) 
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with C^y (p) as in ( |2.13D . Again we can use the same symbol here as for the scalarized form 



factors of the topological charge density since the latter are non-vanishing only for odd 
particle numbers. The above normalization of T^j^ corresponds to the following constraint 
on the scalarized two-particle form factor 

_2Aa2ai 

ha,a, {02, 0l) = JJ. + . . . , 02^01 + lTi. (2.63) 

[02 -01- my 



If we now extend the definitions ( p. 57] ) and ( |2.58| b) for even values of n, we can write 



oo 



P™(^) = Epr^(/^), (2-64) 



k=l 



and define in terms of p^^{fi) as in ( p. 10 ). The various two-point functions can then 



be computed along the lines described in section 2.1. 
Finally the perturbative results for the EM tensor are: 

= i^H'') + 0(a^(A.))} . (2.65) 
In addition, the integral of the spectral density is constrained by ( p.23| ): 

POD 1 

/ d^^p^''{^^) = -, (2.66) 
expressing the fact that the UV central charge of the 0(3) model is equal to 2. 



3. Results for the 0(3) form factors 

In this section we collect our results for various 0(3) form factors, form factor squares 
and their properties and we give an exact expression for the asymptotics of the n-particle 
spectral densities. We refrain from going into the details of the derivation and giving 
proofs, which can be found in 

The most remarkable feature of the NLS model is its integrability. Assuming asymptotic 
freedom one can establish the existence of non-local conserved charges and the absence of 
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particle production ||2^. The matrix elements of the non-local charges between physical 
particle states together with the 0(3) symmetry determine the S-matrix up to a CDD 
ambiguity, and the result agrees with that obtained from the S-matrix bootstrap ||25| . For 
an overview of these issues see also [^]. The 0(3) bootstrap S-matrix has been tested 
(at low energies) in lattice studies |27| and used as an input for the thermodynamic 
Bethe Ansatz (TBA) to compute the exact m/A ratio As a by-product of the 

TEA considerations the CDD ambiguity can be resolved. The construction of the non- 
local charges has been extended to the full quantum monodromy matrix in ||28[. More 
recently the algebraic structure underlying these non-local charges has been considerably 
generalized [EO] and has also been identified directly in the context of the form factor 
bootstrap [^. 

The Form Factor Bootstrap method was initiated in ^ and has been further developed by 
I . Using the exact S-matrix as an input the form factors get constructed as 



Smirnov [31 



solutions of a system of functional equations known as Smirnov axioms. These equations 
recursively relate the n-particle form factors to the {n — 2)-particle form factors and 
(provided the S-matrix also has bound state poles) to the {n — l)-particle form factors. 
The method is reviewed in Smirnov's book |^ and, for the special case of the 0(3) model, 

in m.m. 



3.1 Reduced form factors 

Let us recall the simplest nontrivial case, i.e. the result for the 2-particle scalarized form 
factor of the Noether current 

/M.a.(^^2,^l) = y^(^2-^l)6~, (3.1) 

where 



Here the normalization condition ( p.47|) has been taken into account. The product of the 



2-particle solutions for all possible pairs of rapidity differences. 



^{en,...,e^) = ll^|Ji0,-0,), (3.3) 

i>j 
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will play an important role in what follows. Actually, if there were no internal indices, 
( |3.3| ) would be a solution for the n-particle case. All the following complications (absent 
in integrable models with diagonal S-matrices) are due to the presence of these internal 
quantum numbers. Using (|3.3| ) we can define the reduced form factors Qia^-.-ai i^n, ■ ■ ■ 
for the Spin & Current operators as 

/M„...a. (On, ...,e,) = ^^{9n, ...,e,) gl,^...,^ . . . , ^i) (3.4) 
and similarly for the TC density & EM tensor operators as 

'An -I 

/o,a„-ai {On, ■ ■ • , 6*1) = ^ ^(6'n, • • • , go,a„ . ai (^'n, • • • , ^l) • (3.5) 



The first few reduced form factors are 

9la,a.{02,e,) = e'^'^^'^\ (3.6b) 

90,a2a, {O2, Ol) = , (3.6c) 

02 — Oi — tTT 

^70,a3a,a,(^3,^2,^l) = e~, (3.6d) 



where the normalization conditions ( |2.28| ), ( |2.47| ) and ( |2.63| ) have been used. Note that 



we have not found the normalization condition for the topological charge density operator 
yet. Its lowest reduced form factor (|3.6| d) was chosen here to satisfy the clustering relations 
in the form described below. The still unknown normalization of the physical operator is 
given by ( |2.55| ). 



What makes the reduced form factors useful is the fact that (with the notable exception of 
the 2-particle form factor of the EM tensor (|3.6| c)) they are all polynomial in the rapidities, 
with specified total degree N{n) and partial degree p{n) given by 



N{n) = ^{n^ -3n) + 1 , p{n) = n - 3 + I , 1 = 0,1. (3.7) 

By partial degree we mean the degree in an individual rapidity variable and / = 0, 1 
correspond to the isospin and 1 families, respectively. Expressions for the 0(3) form 
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factors were first written down in M^ . Exploiting their polynomiality and the intertwining 



property, they were recast in p3| , |32| into a form more suitable for practical purposes. 
The fact that suitable reduced form factors are polynomials in the rapidities is a special 
feature of the 0(3) model that fails to hold for the higher 0{N), N > 3 models and 
which considerably facilitates their explicit construction. We have computed the form 
factors for both families of operators (Spin & Current and TC density & EM tensor) 
up to 6 particles. The complexity of these reduced form factors grows very rapidly with 
the particle number, because the number of components (with respect to the internal 
quantum numbers), the number of independent rapidity differences and the degree of the 
polynomials all grow very fast with increasing particle number. For 5- and 6-particles the 
polynomials are already too large to be communicated in print Below we complete the 
list (|3.6|) for up to 4 particles: 

Spin & Current: 



n 



( 9iiO) \ 


( 


^23 \ 


92{e) = 


— 


031 - 2m 


V gsio) J 


V 


ei2 ) 



n 



( am \ 




92{d) 




am 




am 




am 




V aeiO) ) 





+ 5"3a2gaa4ai ^^^^^ ^ ^a.a.^aa.a^ g^^Q^ ^ ^a2a,^aa,as ^^^^^ 



[u — i-njv^v — tTv) 
{u — i7c){u + 2i'K){i'K — v) 
uv{3i7!' — v) 
u{u — m)v 
\ 2m{m — u)v 

( —4:71^ — m{u + v) — {u — vY \ 



+ (^43 - iT^) 



tTT{U + V) — [U — V 

-271^ - Sinv + 
— 47r^ + i7T{u — 2v) — 
2t{^ + i7r(M + 2v) - 2uv 
~i'K{2u + v) + 2uv 
\ -27r2 + i7r(n - 3w) 

where u = 9^2, v = 631. 



+ (^43 - «vr) 







—u 

V — 2i7T 

\ u — V ) 



^For example the 6-particle polynomial of the current has about 3.4 Mbytes, i.e. about 700 A4 pages. 
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EM tensor & TC density: 







{ 


-^21^43 + 27r2 


\ 








^32^41 + ^21^43 - 227r(041 " ^Tt) 




\ am ) 




\ 


-^32^41 + 2z7r(^32 + 


) 



Here we used the shorthand notation 6*^^ for Qi — 6*^. 



3.2 Clustering properties 

A remarkable property of the form factors is asymptotic clustering. If we divide the set of 
rapidities into two disconnected clusters and boost one of them with respect to the other, 
for asymptotically large boosts the form factors factorize into a sum of products of form 
factors corresponding to the two clusters separately. In formulae 

^?,a,-aib„,...6i (ttfc + A, . . . , ai + A, . . . , 

= A'^-^e'^^'^^t.-ai («/^, •••>«!) 9h^-w {Pm, . . . , A) + 0{/\^^-') (3.8) 

and similarly 

5'0,afc-ai6„-6l(afe + A, . . . , «! + A, /^m, . . . , 

= ^'"^-'gla.-a, •••,«!) {Pm, • • • , A) + 0(A^-=^) . (3.9) 



Note that in ( |3.8| ) members of the isospin 1 family are mapped onto themselves, while 
in ( p.9|) members of the isospin 1 family are linked to members of the isospin family. 
Observe also that there is no distinction between even and odd members of the same 
family, their factorization properties are the same. Thus even and odd members of the 
two families are very closely related. We already anticipated this close interrelation (again 
a special 0(3) property) by using the same notation for even and odd form factors and we 
introduced the non-perturbative constants and Aq in (|2.3CI|) and ( |2.55| ) respectively 

in order to have this simple form for ( p.8|) and ( p.9|) . 

For the special case ofA; = l,m = n — 1 the clustering relations read 

9la^-a,{0n, • • • , ^l) = • • • , ^l) + O^-') , 

gO,a^-aA0n, • • • , ^^l) = '^??ra„_,...a. • • • , ^l) + 0(C-^) , (3.10) 
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from which one sees that the reduced form factors are polynomials of partial degree (^ — 2) 
and (n — 3) in the isospin 1 and case, respectively. 

Since the product ( |3.3| ) also factorizes under clustering, the full (scalarized) form factors 
also satisfy clustering relations, which are similar to (|3.8| ) and ( p.9|) . For the / = 1 
family they can be found in Smirnov's book Similar clustering relations were recently 
discussed in . 



The clustering relations closely resemble some classical equations satisfied by our opera- 
tors. For example dividing an even number of particles into two odd clusters, (|3.8|) can 
be interpreted as the quantum counterpart of ( |2.43D , the classical definition of the current 
in terms of the spin operators. (Remember that we are dealing with scalarized objects so 
that all information on the Lorentz structure is lost.) The division of an even number of 
particles into two even clusters, on the other hand, resembles the classical relation 

d,J:-d,Jl = 2gle'^'^j'^Jt. (3.11) 

Finally the clustering of an odd number of particles corresponds to 

d,S'' = -gyS'r^. (3.12) 



Similarly, ( |3.9| ) corresponds to the defining equation (|2.61|) or either of 



T,u = gl{j;Jt - \n,.rjn , q = ^^e^^J^d^S- . (3.13) 



3.3 Form factor squares 

The quantity entering the spectral densities and two-point functions is the absolute square 
of the form factors, summed over the internal symmetry indices. For the reduced form 
factors the corresponding quantities are 

aan---a\ 

Gt\en,...,e{) = ^ E \g,,a^-aAOn,....e,)\\ (3.14) 
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Here are the results for the first few reduced form factor squares: 
G\'\9^) = 1, 
G'i\62,6i) = 2, 

Gf\e,,e,,e,) = 2(61, + ej, + el) + i2n\ 



1 



6-^(^3,^2,^1) = 2, 
6*0^^4, ^3, ^2, ^l) = 2(^43^21 + ^42^31 + ^41^32) 

+A7:\e% + el^ + el, + el, + + el) + 287r^ . (3.15) 

The reduced form factor squares (/ = 0, 1) are symmetric, boost invariant polynomials 
of total degree m? — Zn + 21 and partial degree 2(n — 3 + /). They satisfy the clustering 
relations 

G'p+-)(a, + A,...,«^ + A,/5„,...,A) 

= (/ + !) A^^— ^+2/ _ _ _ ^ G^)(/?^, . . . , A) + 0(A2'='"-^+2^) , (3.16) 



which follow from ( p.8| ) and ( p.9[ ). Again notice that the members of the Spin & Current 
series (/ = 1) are mapped onto themselves under clustering, while the members of the EM 
tensor & TC density series (Z = 0) are linked to the I = 1 series. 

( 3.161 ) is readily translated into the corresponding statement about the full form factor 



squares ( p.31|) and ( |2.57| ) 



+ A,...,ai + A,/3^,...,/3i) 
(12 - Al)n' ^(,)^^^^ ^ ^ ^ ^ . . . , A) + 0(A-^+2.) . (3.17) 



A4-2/ 

Next we discuss two further properties of the polynomials First we present an 

explicit formula for the overall leading terms of these polynomials, i.e. the terms with the 
maximal total degree, which is n'^ — 3n + 21. These overall leading terms are given by 



<^0?oVn, • • • , 6'l) = - I n I Po{(^a{n), • • • , 6'<x{l)) , 

GS5(en,...,^l) = (n^'.l E ^l(^<x{n),...,^.(l)). (3.18) 



22 



Here the summations range over all elements of the permutation group Sn and the func- 
tions Pq and Pi are defined as 

Po{9n, . . . , 9l) = , Pi(6'„, . . . ,6*1) = -2 /j2 n2 ' (^-l^) 

respectively. These overall leading terms G^^q by themselves satisfy the clustering relations 
( |3.16| ). An equivalent expression for the leading terms ( |3.18| ) was obtained by H. Lehmann 
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The second property concerns the analytic continuation of the squares G["^ to complex 
rapidities. For real, physical rapidities the polynomials G^'^^ are real-valued and positive. 
Being polynomials one can also evaluate them for complex rapidities 6i. In particular 
they turn out to vanish if two consecutive rapidity differences are both equal to ivr 

G\''\9n-2 + 2m,9n-2 + i7^,9n-2,...,ei) =0, n>4-/. (3.20) 

G^^^ vanishes also at all other points in rapidity space obtained from the one in ( |3.20| ) by 
permutation of the arguments. 

Since the form factor squares - not the form factors themselves - are the objects enter- 
ing the physically interesting quantities through the spectral representation it would be 
desirable to obtain them directly, without having to go through the tedious procedure of 
first computing the form factors and then squaring them. Indeed, for the first few n the 
information collected in this subsection is sufficient to determine G^""^ directly. However, 
for particle numbers 5 and 6, although these constraints - permutation symmetry, boost 
invariance, given total and partial degrees, given overall leading terms, clustering and the 
vanishing at points (|3.20|) - almost uniquely determine the solution, but not quite. So in 



these cases we had to work out the solution the tedious way. We hope to find additional 
general properties of the squares that will determine them completely. In the appendix we 
list the results for the form factor squares for the EM tensor & TC density and Current 
& Spin series up to 6 particles. 



3.4 Asymptotics of spectral densities 

Let us now turn to the computation of the spectral densities for the four operators under 



23 



consideration. The n-particle contributions are 

(n)/ X _ dUi... dUn-l ^(n), 



pr\f^) = I \;^:^y:.r' Fi-\u) ^(/.-m(")(^)), n>2, 



Ft\u) = -^{3-2l)Gl-\u)\^{u)\\ (3.21) 

The integration in ( |3.21| ) can be done numerically. The graph of the resulting n-particle 
spectral densities is roughly 'bell-shaped': Starting from zero at /i = mn they are strictly 
increasing, reach a single maximum and then decrease monotonically for large /i. Some 



typical plots can be found in 



Of particular interest is the fi ^ oo asymptotics of the spectral densities. One can show 
that an important consequence of the clustering relations ( |3.17|) is the following expression 
for the asymptotic form of the n-particle spectral densities 

where the constants A^"^ are computable from the integrals of the lower particle spectral 
densities. The integrals are 

cS") = dfipt\fi) = I (0'"/o" • • • d''r^-lG^^\u) \^{u)\' , n>2, (3.23) 
and in terms of them the decay constants A\"'^ are given by 

Q n—1 

^(n) ^ « ^(n) ^ ^ ^ ^{k) ^in-k) ^ ^^^^ ^3^24) 



k=l 



if ( |3.23|) is supplemented by the definition := ir/A. Note that (|3.24|) implies that the 



a\"'^ are strictly increasing with n. This implies that for sufficiently high energies the n- 
particle contribution overtakes the {n — 2)-particle contribution, i.e. p*^"'^(/i) > p'-"'"^^ (/i) , 
for /i > where p'^^\iin) = P*-"'~^H/^n); although the maximum of p^"^(/i) is expected to 
be much smaller than the maximum of p*^""^-* (/i) (c.f. section 4). This is the "crossover 
phenomenon" observed in for low particle numbers. Although ( |3.22| ) is an exact 



asymptotic equation, the way how the asymptotic behavior is approached is highly non- 
uniform in the particle number. That is to say, with increasing n one has to go to larger 
and larger p in order to make the right hand side of ( p.22| ) a good description of the 
function p["^(p). 
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A consequence of the crossover is that the asymptotic behavior of the exact spectral 
density (being the sum of all even/odd n-particle contributions) cannot be obtained by 
naively summing up the right hand sides of ( p.22|) , which in fact would be divergent. This 
divergence signals that the infinite series has a more singular high energy behavior than 
each of its terms separately. Indeed, according to PT, the high energy behavior of the full 
spectral densities is 

O = EM & top: p(/i) 

/i 

O = spin & curr: p(/x) ~ 

where the values of the various decay constants can be read off from the equations ( |2.40| ), 
( |2.51| ), ( |2.60| ) and ( p.65| ). Clearly a finite number of terms, each decaying according to 
( |3.22| ) can never produce the more singular UV behavior of ( ^.25] ). It is therefore a major 
challenge to develop re-summation techniques for the spectral resolution and, for the 
reasons outlined in the introduction, compute the UV properties independent of PT. We 
shall meet this challenge in section 5. 



/ In In p \ 



1 + 



In fi 



(3.25) 



4. Results for the two-point functions 

Knowing the form factor squares the evaluation of the n-particle contributions to the 
spectral densities and the two-point functions is in principle straightforward. For the 
spectral densities the integrations in ( |2.32| b) and ( |2.58| b) can be done numerically to good 



accuracy. Throughout we used an accuracy of 10^^ for all numerical computations. For 
comparison with PT and MC data it is useful to consider the Fourier transform I{p) of 
the two-point function, which can be computed from the spectral density via (|2.10| ). Of 
course I{p) again decomposes into a sum of n-particle contributions of which in practice 
only the first few are known explicitly. It is important that an intrinsic estimate for the 
error induced by this truncation can be made, i.e. one which does not rely on comparison 
with other techniques. 

To discuss this qualitative error estimate let us first note some general features of the 
n-particle contributions to the spectral densities. As remarked before the graph of an n- 
particle spectral density is roughly bell-shaped: Starting from zero at /x = mn it is strictly 
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increasing, reaches a single maximum and then decreases monotonically, the /i ^ oo 
asymptotics being given by ( |3.22| ). With increasing particle number n the values of 
the maxima rapidly decrease. Generally speaking the maximum of p''"^(yu) is smaller 
by 1.5 to 2.5 orders of magnitudes compared to the maximum of p*^"~^-'(/i), while its 
position is shifted to higher energies. Nevertheless at sufficiently high energies the n- 
particle contribution overtakes the {n — 2)-particle contribution, i.e. p^"'\fi) > p''"^^-*(/i), 
iff /i > fin, where p^^\fj.n) = P^^~'^Kl^n) For n < 6 the positions and the values 

of the maxima are listed in Table 1 below. The results for the points of intersection 
{(Xn/'nT', p^'^~'^\iin) = P^^K^n)) are as follows: 



-6^ 

5 



Spin: n = 5: (1.0 ■ 10^ 4.9 ■ 10" 

Current: n = 4: (1.6 ■ 10^, 3.9 ■ 10"^) , n = 6: (1.0 ■ 10^ 3.9 ■ 10"^) , 

TC: n = 5: (7.0 ■ 10^, 3.8 ■ lO^^) , 

EM: n = 4: (1.6 ■ 10^, 4.7 ■ lO'^) , n = 6: (4.6 ■ 10^ 2.5 ■ 10"^°) . (4.1) 

It is natural to assume that the general trend depicted by these numbers continues to 
hold at higher particle numbers; a quantitative form of this assumption will be presented 
in section 5. Here already the qualitative features (/i„ increasing; maxima decreasing) 
are sufficient to conclude that the point of intersection /i„ provides an intrinsic measure 
for the quality of the approximation made by truncating the form factor series at the 
n-particle term: Since Pn+2 ^ Pn the (n + 2)-particle contribution can safely be ignored 
up to energies p^Pn in the sense that the correction to p^'^\p) + . . . + p^^\p) for p^pn 
should be at most a few per mille. Thus, the form factor series truncated at 6 particles 
should provide results accurate to within a few per mille for the spectral density at least 
up to energies 7 ■ lO'^m for Spin and TC and 5 ■ lO^m for Current and EM. It is not hard 
to see that for the functions I{p) related to p(/i) by the Stieltjes transform ( |2.10| ) this 
implies accurate results for somewhat smaller ranges, about p^lO^ for the Spin and TC 
cases and about p^lO^ for the Current and EM tensor. Both of these energy /momentum 
ranges exceed that accessible through MC simulations by several orders of magnitudes. 

Below we shall restrict attention to the Spin field and the Noether current. In both cases 
it is convenient to consider the low- and the high energy region separately. In the low 
energy region non-perturbative effects are expected to become important. Here Monte 
Carlo simulations provide an alternative non-perturbative technique to probe the system 
27| , |35| , 0, |l^]. Simulations for the two-point functions were made [0 using a Wolff-type 



cluster algorithm p2| on a 460 square lattice at inverse coupling (3 = 1.80 (correlation 
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length ^ = 65.05). A comparison between form factor results, MC data and PT in this 
low energy region is shown in Figures |I| and ^ for the Spin and Current case, respectively. 



2 2 
P I[p ] 




Figure 1: Low energy region of spin two-point function. Comparison: Form factor ap- 
proach, Monte Carlo data and 2-loop perturbation theory; p^J(p^) plotted against p/m. 
The normalization of the PT curve is fixed according to ( [^.41| ). 

The agreement between the MC data and the form factor results is excellent. In the Spin 
case the statistical errors of the MC data are less than the size of the dots in Figure 
|l[ Above 30m they have a slight tendency to lie below the form factor curve, which 
(comparing e.g. with the data at (3 = 1.7) can be attributed to the still finite correlation 
length. The 1 -|- 3 ff curve is left out in Figure 1 as it coincides with the 1 -|- 3 -|- 5 ff 
curve below 30m and (incidentally) with the PT curve above 40m. In the Current case 
the statistical errors are larger but within the errors the agreement with the form factor 
curve is perfect. One also sees that for energies between 30 and 45 the PT curve runs 
almost parallel to the MC data in Figure Without the guidance of the form factor 
result one would thus be tempted to match both curves by tuning the Lambda parameter 
appropriately. Doing this however, the Lambda parameter comes out wrong by about 
10% (from below), as compared with the known exact result |1^. Generally speaking one 
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sees that a determination of the Lambda parameter (accurate to within 1% say) from MC 
data and (2-loop) PT alone is difficult because some assumption about the onset of the 
2-loop perturbative regime enters. It remains difficult even after 3-loop effects are taken 
into account To investigate the onset of a perturbative regime let us now consider 
the high energy regime of the same quantities. 

In Figures ^ and ^ again the Fourier transform /(p^) of the 2-point function of the Spin 
field and the Noether current are shown, computed once in 2-loop PT and once via (|2.10|) 
by truncation of the form factor series. In the spin case the normalization of the PT 
curve is fixed according to ( p.41| ), so that as in the current case no free parameter enters 
the comparison. Let us define the perturbative regime to be the energy range for which 
the 2-loop PT and the truncated form factor curve coincide within 1%. In both the 
Spin and the Current case a large perturbative regime is found to exist. Having fixed 
m/A and Ai no adjustable parameters enter the PT results. The very existence of such a 
perturbative regime thus supports the proposed value of m/A as well as the exact value of 
the normalization constant ( |2.41|) . In the Spin case the 2-loop perturbative regime is about 
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Figure 3: Spin two-point function. Comparison: Form factor approach versus 2-loop 
perturbation theory; logplot of against p/m. The normahzation of the PT curve 

is fixed according to (|2.41| ). 

50<,/i/m<,5000 and in the current case it is about 100<,/i/m<,2 ■ 10"^. The lower bound 
of this interval will remain unaffected when higher particle form factor contributions are 
taken into account and thus is a genuine feature of the 0(3) model. Notice however that 
the onset of the 2-loop perturbative regime occurs at much higher energies (50 — 100 
times the mass gap) than is sometimes pretended in the 4-dimensional counterpart of this 
situation. Probably this should be taken as a warning in 4-dimensional gauge theories. 

The upper bounds /x ~ 5000 for the (1 + 3 + 5)-particle Spin curves and /i ~ 10"^ for 
the (2 + 4 + 6)-particle Current curve are fully consistent with the intrinsic estimate for 
the validity of the truncated form factor results given at the beginning of this section. 
These upper bounds of the perturbative regime, of course, are expected to move up as 
higher particle form factor contributions are taken into account. Indeed, the hypothesis of 
asymptotic freedom is that they approach infinity, so that the entire high energy regime 
is perturbative. As emphasized in the introduction, the issue of asymptotic freedom can 
only be addressed if the high energy properties of the theory can be explored by means 
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Figure 4: Current two-point function. Comparison: Form factor approach versus 2-loop 
perturbation theory; logplot of p^I{p'^) against p/m. 

of a rehable technique independent of PT. Guided by the form factor results we formulate 
in the next section a conjecture about the structure of the n-particle spectral densities 
for general n. Assuming the validity of this conjecture the extreme UV behavior of the 
two-point functions can be computed. 



30 



5. Scaling Hypothesis for Spectral Densities 

The 77,-particle contributions to the spectral densities appear to follow a remarkable 'self- 
similarity' pattern. Suppose the scalarized n-particle spectral density of one of the four 
local operators considered to be given. Can one - without actually doing the computation 
(which for large k is practically impossible) - anticipate the structure of some k > n 
particle spectral density? We propose that the answer is affirmative and for large n and 
A > 1 is given by 

P^'^-Hp) - ^ ipM'/''"" p(-) (m(p/m)VA-") , (5.1) 

where m is the mass gap and the exponents 7 and a are given below. Taking A = for 
example yields a candidate for the {n + 2)-particle spectral density. In the following we 
shall first give a precise formulation of the scaling law ( p.l| ) and present evidence for it. In 
section 5.2 we then promote (|5.1|) to a working hypothesis and explore its consequences. 
The results obtained have an interesting interplay with both PT and non-perturbative 
MC data. All pieces of information are found to be consistent, which we interpret as 
further supporting the hypothesis. 

5.1 Evidence for Self-similarity and Scaling 

Since the leading asymptotics of the spectral densities is given by it is useful to 
consider fxp^^'^fj.) instead of p^"-'(yu). Further a logarithmic scale is convenient so that we 
are lead to introduce 

i?f"^(x) := me^p{"^(me^) , lnn<x<oo, / = 0, 1 . (5.2) 

Here Z = 0, 1 as before correspond to the EM tensor & TC density and Spin & Current 
series, respectively. The graphs of these functions are again roughly 'bell-shaped': Starting 
from zero at x = Inn they are strictly increasing, reach a single maximum at some 
X = ^l"^ > Inn and then decrease monotonically for all x > ^l^K The position ^l^^ of 
the maximum and its value m/"^ := i?["^(^|^"^) are two important characteristics of the 
function, and hence of the spectral density. For n < 6 these data are collected in Table 1. 

The content of the scaling law ( |5.1| ) is most transparent if the ordinate and abscissa of 
the graph are rescaled such that both the value and the position of the maximum are 
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Table 1: Positions and values of the maxima 



n 


An) 
Co 


1 n4 ?i /r(n) 
10 Mg 


^1 


10 Ml ' 


2 


0.953 


317.3 


1.155 


27.33 


3 


2.726 


27.95 


3.613 


10.76 


4 


4.720 


7.870 


6.631 


6.736 


5 


6.945 


3.134 


10.06 


4.911 


6 


9.344 


1.517 


13.79 


3.867 



normalized to unity. Thus define 



(n) 



(lnn)/^{"^ < ;z < oo , / = 0, 1 . 



(5.3) 



In order to have a common domain of definition we set yI'''\z) = for < ;z < (lnn)/^J"\ 
The proposed behavior of the spectral densities is as follows: 

Scaling Hypothesis: 



(a) (Self- similarity) The functions Yi^^\z), n > 2 converge pointwise to a bounded 
function Yi{z). The sequence of k-th moments converges to the k-th moments of 
Yi{z) for k + 1 = 0, 1, i.e. 



\im Yt\z) = Y,{z) , z>0 



lim / dzz''Yf''\z) 



dzz''Yi{z) . 



(5.4) 



(b) (Asymptotic scaling) The parameters C.i'^^ and M^^'' scale asymptotically according 
to powers of n, i.e. 



(5.5) 



Note that the convergence in (a) is weaker than uniform convergence. Feature (a) in 
particular means that for sufficiently large n the graphs of two subsequent members 
y/"-'^(z) and ^/"^(z) should become practically indistinguishable. This appears to be 
satisfied remarkably well even for small n = 4, 5, 6, as is illustrated in Figures ^ and 
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Figure 5: Illustration of the self-similarity property of the rescaled / = spectral densities. 
The plots show Y^"'\z) (dashed) compared with Yq'^~^^\z) (solid) for n — 2,3, 4, 5. 



To substantiate proposal (b) let us prepare some further characteristics of the functions 



Y[z] 
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Figure 6: Illustration of the self-similarity property of the rescaled / = 1 spectral densities. 
The plots show y}"'\z) (dashed) compared with Yi^~^^\z) (solid) for n = 2, 3, 4, 5. 
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rI"'\x), / = 0, 1: Its integral, its first, respectively minus first moment, and the strength 
of the asymptotic decay. In formulae 



127r 



dxR 



X) 



(n) 



R^-\x) 



in), 



X] 



A. 



(n) 



X 



oo 



x^ 



(n) 



dxRi 
dx 



X] 



r^-^Rt^ 

Jo X 



X) 



(n) 



X 



X^ 



oo 



(5.6) 



In practice of course one will evaluate the constants c["'' and h\"'^ directly as {n — l)-fold 
unconstrained integrals, as in (|3.23|) , rather than first computing the spectral densities 
and then evaluating their moments. For n < Q the results for these constants are listed 
in Table 2. 

Table 2: Integrals, ± first moments and decay constants 



n 


(n) 
Co 


102^") 


An) 
^0 


An) 
Cl 


lO/i^) 


^(n) 


2 


1.6027 


8.231 


0.727 


1.009 


3.958 


1.938 


3 


0.4104 


5.251 


1.866 


1.140 


1.693 


4.977 


4 


0.1943 


4.140 


3.308 


1.242 


1.054 


8.821 


5 


0.1117 


3.430 


5.006 


1.327 


0.761 


13.35 


6 


0.07185 


2.930 


6.938 


1.400 


0.594 


18.50 



In addition one has = tc/A. The decay constants A^""^ are difficult to compute directly 
since the truly asymptotic behavior of the spectral densities sets in only at astronomically 
large energies /i ~ me^°"^''^. Using the exact relations p.24| ) they can however accurately 
be computed from the integrals Ci^\ 



Making use of the convergence ( |5.4| ) one can deduce scaling laws for the quantities ( |5.6|) . 
To this end we introduce the corresponding quantities for the universal shape functions 

PCO POO 

c*. = 127r / dzYo(z) , c? = / dzYi(z) , 

JO Jo 
rco rco fly 

h*= dzzYoiz), h\= —Yr{z), (5.7) 

Jo Jo z 

A* A* 

Yo{z) J , z ^oo , Yi{z) ^ , oo . 

z^ z^ 
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For large enough n we can approximate b!(^\x) pointwise and w.r.t. the moments by 

rill ' Uirii+°' 



Combining ( |5.5|) - ( |5.8D one finds 

4"^ ~ /iS^O^O ^2+2^0-70 ^ /^W _ /^*]Vf^ ^-71 ^ (5.9) 



-71 



Two important constraints arise from the non-hnear relations ( p.24| ) for the asymptotic 
decay constants Af^\ Combined with ( |5.9| ) this yields the following relations among the 
parameters 

70 = 3 + 4ao - 2ai , A^Mo^o = ^^IMi^? • (5.10b) 



(|Oqb) follows trivially from ^S"^ = To arrive at ( p).l(J| a) we approximated 



n-l 1 

Y^[k{n - k)f ^ n^^+' / 
k=i -^0 



which is valid for large n. The exponents 7; are determined in terms of ao and ai, which 
will later turn out coincide 

aQ = ai =: a . (5-11) 



Re-inserting ( p.lO| ) and ( |5.11| ) into ( [5.9| ) one arrives at a simplified set of scaling relations. 
In particular 

^ n , hi:= hlMi , 



ht"^ hi 
1 1 



/iS") ho 



n , ho:= hlMoi^ 



2 

) 



[4")]-V(2+.) ^ ^ (c*Moeo)-^/(^+") , (5.12) 

are predicted to scale linearly with n. The same holds for the following two ratios, which 
provide a convenient way to determine the remaining exponent a from the slope of the 
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linear function 

—r-\ — = 7—\ — ~ TT— n . 5.13 

3 [c(")]2 r(2a + 2) ^ ' 

The predicted linear scaling in ( ^.121) , ( ^.13 ) appears to be satisfied fairly well even for 
small n = 4, 5, 6 as is illustrated in Figure |^ 




Figure 7: Two-parameter linear fits for the n-dependence of various parameter combina- 
tions. 

The fits are always two-parameter linear fits (slope and ordinate) on the largest three 
data points n = 4,5, 6. Using (|5.13| ) the exponent a comes out to be 



a ^ 0.273 . (5.14) 

Using this value in the fitting of [cq"'']"^/''^"''"-* the predicted linear scaling is confirmed. 
The quality of the l/ho^^ fit confirms ( p. 11 ). For later use let us also note the values of 



the slopes in (|5.12| ) 

l//ii = 3.674, l//io = 4.988, Cq = 0.5643 . (5.15) 

Similar fits can be made for the values and the positions of the maxima, given numerically 
in Table 1. They are similarly convincing; we refrain from displaying them. 
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5.2 Summation of n-particle contributions 



Having at hand candidate expressions for tlie n-particle spectral densities as given in ( ^.8|) 
one can evaluate their sum. Suppose that for n < uq the n-particle spectral densities have 
been computed explicitly; in our case no = 7. Decomposing the sum 

riQ — 1 oo 

Riix)= E E (5-16) 

n=2-l n=no 

we wish to evaluate the second term. Using ( ^.8|) the sum can be expressed in terms of a 
suitable integral transform of the shape function Yi{z). Here we shall need only the large 
X asymptotics of the sum, in which case it is sufficient to approximate the sum over n by 
an integral. One finds for the large x behavior 

where the corrections are at least one power down in x. Clearly the physical spectral 
densities can be evaluated similarly; summing only over odd/even particle numbers will 
simply result in an extra factor 1/2 in the asymptotic expression ( [5.17| ). In particular 
in the isospin 1 case the exponent 71 = 1 yields the asymptotics Ri{x) Ai. One can 
show that the asymptotics of Ro{x) must always be down by two orders in x, compared 
to that of Ri{x). This is related to, although not a direct consequence of ( |3.22|) . One 
concludes that the x ^ 00 asymptotics in the isospin zero case must be Ro{x) Aq/x"^. 
Comparing with ( |5.17D this requires ao = ai, as claimed before. Moreover one obtains 
the relations 

^EM _4^top _ -A_^ , (5.18b) 

which have various interesting consequences. ( p.l8| a) relates A^^^^ to A^"^™. This is in- 
teresting because A'^^^^ is accessible to perturbation theory, while A^^^^ is some previ- 
ously unknown non-perturbative constant, which via ( p.38|) , ( p.40|) determines the short- 



distance normalization of the Spin two-point function. On the other hand both constants 
also get expressed in terms of the scaling hypothesis. Using ( p.l4| ) and ( |5.15| ) one finds 



{A^'''')sH = 1/9.3553 = 1.007/37r. Since perturbation theory gives (A="")pr = l/37r 
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one sees that not only the functional form of the asymptotic behavior obtained from the 
resummation of the spectral resolution agrees with that predicted by PT, but also the 
numerical value of the leading coefficient agrees with an accuracy better than 1%. De- 
pending on the viewpoint one can interpret this as supporting the validity of PT and/or 
our scaling hypothesis [|. Accepting the PT value for A^^^^ we get the exact prediction for 
the non-perturbative constant Ai = A'^p™ anticipated in ( |2.41|) 



Ai = A'^P''^ = — = 0.135095 . (5.19) 
37r"' 



On the other hand A'^p™ has been evaluated by Monte Carlo simulations [|I^ with the 
result (y4^P™)A/c = 0.135(2), convincingly close to the proposed exact result. The similar 
behavior of the Spin and Current two-point functions and spectral densities for large 
energy is again an 0(3) speciality, as can be seen from ( |2.42| ). 

A similar pattern underlies ( |5.18| b). Now both ^4^^ and y4*°P can be computed in PT. The 
results are given in ( 2.60| ), (|2.65|) and yield in particular {A^^)pt = 4:{A^"^)px- The non- 



perturbative result now is that this transfers to relations for the exact spectral densities 
at all energy scales 

OO 1 oo 

k=l ^ k=l 

which is the announced result ( |2.56D . Consistency then requires that the numerical value 
for A^^'^ computed via the scaling hypothesis coincides with the perturbative value. Com- 
bining (|5.14|) , (|5.15|) and ( |2.65|) one finds 



(A^^i)^^ = 1/12.700 = 0.989/47r , (A^''')pt = l/47r . (5.21) 

Note that here the matching between full and perturbative djTiamics already concerns 

a 1-loop coefficient and is accurate to 1%. Again this cuts both ways, supporting the 

validity of PT and/or our scaling hypothesis and thus also the proposed normalization 

constant Aq = 1/4 of the topological charge density operator. 

"'^The extrapolation based on the scaling hypothesis contains a (neglegible) numerical error as well as 
a systematic error. The latter is due to the uncertainty introduced by doing the fits for moderately large 
particle numbers n=4,5,6. In principle there may also be subleading (not powerlike) terms in the scaling 



laws (5^). A rough feeling for the size of the systematic errors can be obtained by ad-hoc changing a in 



( pl| ) to a = 0.25. One finds (^™")sif = 1.026/37r; {A^^)sh = 1.008/47r; csh = 1.998. 
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Finally consider the central charge as defined in ( |2.18| ). In the spectral resolution it again 
decomposes into a sum of n-particle contribution. The first few terms are given in Table 2, 
i.e. 



ci'^ = 1.603 , c^o^ = 0.194 , = 0.072 , cj,"^ = 127r H d^i p^^\^i) . (5.22) 

JO 

The remainder of the series can be evaluated by means of the scaling hypothesis. Using 
(|57[5| ), ( CT ) one finds 



= cS,'^ + cS,'^ + 4'^ + Y.{2kco - 0.200)~(2+") = 1.997 . (5.23) 

fc=4 

Here we also used the ordinate of the linear fit. (Tree level) PT predicts cpt = 2, cor- 
responding to the two unconstrained bosonic degrees of freedom. The form factor com- 
putation shows that this is compatible with the non-perturbative low energy dynamics 
of the model and provides further support for the scaling hypothesis (see the footnote on 
page 38). An alternative non-perturbative consistency check is provided by the thermo- 



djTiamic Bethe Ansatz |38l] and also yields c = 2. Let us point out that the value c = 2 
does not carry much information about the nature of the UV limiting CFT. For a number 
of reasons it cannot be the Gaussian c = 2 model. 



5.3 Towards the exact two-point functions 

The scaling hypothesis also allows one to numerically compute the two-point functions 
at all energy/length scales. In coordinate space this simply amounts to performing the 
integral in the spectral representation, once the full spectral density is available via ( |5.16| ) 
and ( ^.8|) . In momentum space one can combine eq. ( p.lO|) and (|5.8| ) to evaluate the 
sum over all n-particle contributions to the functions I{p), basically giving the Fourier 
transform of the two-point function. The results will be described in p9 . 



6. Conclusions 

In this paper we studied the two-point functions of the four physically most interesting 
operators in the 0(3) NLS model. We have chosen the 0(3) model because it is the 
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simplest 2-dimensional theory that resembles in many aspects QCD. Concerning the on- 
shell features various exact results, in particular the exact two-particle S-matrix are known 



25, 24], due to its remarkable integrability properties. On the other hand there exist a 



large number of MC studies of the model \U7i |22| , pq , |Tg, not making use of this 
integrability. Here we employed the form factor bootstrap method for integrable QFTs 
to compute off-shell features that can be compared both with PT and with MC data. 
A simplifying feature of the 0(3) model (as compared, say, to the higher O(A^) models) 
is that the computation of the form factors is relatively straightforward due to their 
essentially polynomial nature. Once the form factors are known, the two-point functions 
can be evaluated in terms of a spectral resolution; that is as an infinite sum where multiple 
integrals of the squares of the n-particle form factors enter, the sum running over all 
particle numbers n. In practice one has to truncate this infinite sum after the first 
few terms. In models with a diagonal S-matrix and a power-like approach of the UV 
limit, application of this technique showed extremely fast convergence of the truncated 
sum |jlO|, |Tl], [l^.Q In an asymptotically free theory one expects less fast convergence, 
because the approach to the asymptotic form is only logarithmic. Thus, in order to really 
check the validity of PT (which has been questioned [0, |[) we needed to extrapolate 
our results beyond the largest particle number, presently 6, for which we carried out the 
computations explicitly. The extrapolation was based on a novel scaling hypothesis. Using 
this hypothesis we have been able to compute the extreme UV properties independent of 
PT. Remarkably we found very good agreement with PT: Within 1% for a RG improved 
tree level coefficient in the case of the Current and also within 1% for a RG improved 
1-loop coefficient in the case of the EM tensor. This is very strong evidence for PT. The 
agreement between our extrapolated data and PT also indirectly confirms the proposed 



exact value of the perturbative A parameter [|I^. On the basis of the evidence presented 
in section 5 we regard this scaling hypothesis as very plausible, although an analytical 
proof remains to be found and would be highly desirable. On the other hand, the nature 
of this hypothesis is completely independent of the assumptions underlying the derivation 
of the usual RG improved expressions based on PT, so that the agreement of the results 
is still remarkable. 

A further result is the exact determination of two, previously unknown, non-perturbative 
constants Aq and Ai. Using the values of these numbers, we can write the exact 3-particle 



^Nevertheless we expect the convergence stih to be non- uniform in the sense that the UV behavior of 
the infinite sum is different from that of each partial sum. 
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matrix element of the properly normalized topological charge density operator q{x) as 



7 

7r2 r 



(0| g(0) las, 9s; a^, ^2; a,, 9,) = — [m^ - M^^) (^3, 02, 9^f\ ^{9^, ^2, ^1) e'^^'^^'^^ , (6.1) 

where M^^\9-i,92,9i) is the three-particle invariant mass. Further the short distance 
expansion of the Spin two-point function is now unambiguously fixed to be 



1 2 

(In mr) -|- O (In In mr ■ In mr) 



37r3 



(6.2) 



where m is the mass gap and r = Jxf + X2- Note that the values of the overall constants 



in (|6.1| ) and ( |6.2|) follow from the scaling hypothesis alone, we do not need here the results 
of any numerical fitting procedure. 

Finally we remark that our scaling hypothesis may be viewed as a 2-dimensional analogue 



of the KNO scaling |^ in QCD. We also expect that scaling hypotheses of a similar type 
can be formulated in many other 2-dimensional models, giving simultaneous access to 
their off-shell properties at all energy/length scales. 
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Appendix: List of form factor squares 

Here we list the results for the form factor squares for the EM tensor & TC density 
and Current Sz Spin series up to 6 particles. The corresponding Mathematica files can 
be obtained from the authors upon request. The squares are boost invariant symmetric 
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polynomials in the rapidities and therefore conveniently described in terms of a basis in 
this space of polynomials. 

Let P^''\N,p) denote the space of homogeneous symmetric polynomials in n variables 
that are of total degree and partial degree p. By partial degree we mean the power 
with which an individual variable enters. A basis for P'^"-\N, p) can be obtained as follows. 
For fixed n let (T]^\ . . . , cr^'*-' denote the elementary symmetric polynomials in 6i, . . . ,9n, 
i.e. 

E (A.3) 



a 



(n) 
k 



ti<...<tk 

Let A = (Ai, . . . , Ap), Ai < A2 < . . . < Ap be a partition of N into p parts less or equal 
to n, i.e. J2i\ = 1 < Aj < n, I < i < p. Running through all those partitions, the 
assignment 

(\ \ \ , Jn) (n) 

provides a basis of P^'^\N,p). However these functions arc not boost invariant and one 
would like to have a description of the boost invariant subspacc P/"J(A^, p) of P^'^\N,p). 

(n) 

Boost invariance can be incorporated by switching to symmetric polynomials defined 
by 

ri''\e,,...,en)^a^^\e,,...,en), 2<k<n, 

rj"^ = -4'^ , e, = e, - -{Oi + ... + 0n). (A.4) 

n n 

All monomials in the r;["^'s (k > 2) are manifestly boost invariant. The price to pay is that 
the partial degree is no longer manifest. For a monomial in the cr^^''s the partial degree 
is simply given by the number of factors. This is no longer the case for monomials 
in the T;["^'s. In fact, inverting the relation 

= + E : r,, k^l,...,n (A.5) 

(here and henceforth the superscripts (n) are suppressed) one sees that = (1~ 
k)n^^ + . . . + (jfc. Thus for a monomial in the r^'s the total and the partial degree coincide. 
Linear combinations of such monomials however may have a partial degree that is less 
than their total degree. 

The general structure of the polynomials G^"'* {9) can now be described as follows. 
Gt^ (6) = (6) + Tracts (6) + ...+ 7r^Gg(^) , e p£ (N - 2k, p) . (A.6) 
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The overall leading terms gI"q{6) are given in ( |3.18| ). Below we list the results for the 
polynomials Gl^\d), n < 6 for the Spin & Current and EM tensor & TC density series: 

Spin & Current 

= 1 , Gf\9) = 2 , Gf\9) = 12[-T, + n'] , 

G'i\e) = -4[6r| + 9t^ + AOT2T4] + 8tt'^[25tI + Uu] - 44871^2 + 2727r6 , 

Gf\e) = [-48r|r| - 72r| + 144x2^4 + 52x3X3^4 + 352r^T^ - 640r| + 
+660x2^3X5 + 3200x3X4X5 + 4500x2x2] 

-47r2[36x2^ - 64x2^1 + 692x|x4 + 295x|x4 + 80x2X4^ + 2925x2X3X5 + 4375x5^] 

+ 167r^ [160x2^ + 97x2x| + 909x^x4 - 20x2 + I95OX3X5] 

-167r^[978x| + 380x2 ^ 1965^2^4] 

+647r^ [678x2 ^ gg^^j _ 55^20 vr^S + 24960 ti^^ , 



.3 



Gf\e) = A[~2AtItItI - 36X3V42 + 72X2VI + 26x2X32x1 + 176x22x4^ - 320x, 
+72x2^1x5 + 108x|x5 - 228X2V3X4X5 - 57x2x^x4X5 - 446x22x3X42x5 + 2208x3x|x5 
+ IO8X2V52 - 801x22x32x52 + II82XIT4X52 - 3615x32x4X52 - 430x2X42x52 + I8OOX2X3 
+5625x5^ - I8OX2V32X6 - 62IX2X3V6 + 480x|x4X6 + 2082x22x32x4X5 - 2344x^X42x6 
-1710x32x42x6 - 128x2X4^6 + 3204x|x3X5X6 + 3753x|x5X6 - 744x2X3X4X5X6 
+ 14490x22x52x6 - 34950x4X52x6 - 288X2V2 + 13230x2X32x52 - 18696x22x4X62 
+40464x42x02 + 48870x3X5X02 + 61992x2X0^] 

-47r2 [72x^x1 + 108x36 - 384X2V32X4 - 4O8X2X3V4 + 6OOX2V42 + 2QtItItI 
+2248x2^4^ + 1846x32x| - 1504x2X4^ - 264x3^3X5 - 582x22x^X5 - 3518x^x3X4X5 
-6525x|x4X5 + 8258x2X3X42x5 + 5574X2V52 - 11685x2X32x52 + 22220x22x4X52 
-14450x42x52 + 32625x3x| + 1056x^X6 + 41 10x^X32x6 + 89IX3V6 + 2768x3^x4X5 
+23766x2X32x4X6 - 5U96t^t^tq + 30848x|x6 + 86I9OX2V3X5X6 - 54090x3X4X5X5 
+78600x2X52x5 - 47688x2Vg2 + 136404x32x02 + 51432x2X4X52 + 378000xg3] 

+ 167r^[-126x|x32 - 24x22x3^ + 336X2V4 - 1272x|x32x4 - IO2OX3V4 + 38I8X2V42 
+3769x2X32x42 + 426x32x4^ + 376x4^ - 6I2X2V3X5 - 6291x2x|x5 + 6302x22x3X4X5 
+876x3X42x5 + 21468X2V52 - 135x32x52 + 23185x2X4X52 + 5640x|x6 + 29352x22x32x6 
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-29108x2^4X6 + 19149x3^4X6 - 41452X2X4V6 + 142074x2X3X5X6 + 33075X5V6 
-93966x2V| + 159882x1X6^] 

-87r^[408Tj - 4482x2^1 - 6564x2X3^ + 18276x|x4 + IOOI6X2V3V4 + 45484x|x| 
+203OIX3VI - 5116x2x| + 24628x|x3X5 - 37755x|x5 + 117453x3X3X4X5 + 264035x1x5^ 
+72675x4X5^ + 69236X2V6 + 381147x2x|x6 - 507872x|x4X6 - 17628x|x6 

+623655x3X5X6 - 418404x2x|] 

+327r^ [2964x2^ - 3444x1x2 - 6138x3^ + 409OOX2V4 + 30984x2x1x4 + 26211x2V| 
-1048x| + 69483x2^3X5 + 62973x3X4X5 + 181365x2x| + 33996x|x6 
+182187x3^6 - 353994x2X4X6 - 129717x6^] 

-327r^°[31348x| + U205r^ri + 177930x^X4 + 60772x^X4 
+18686x2x| + 237656x2x3x5 + 184975x5^ - 6386x^x6 - 357082x4X6] 

+327r^2 [168412x2^ + 78879x2x| + 4I8856X2V4 - 8772x4^ + 26O8IIX3X5 - 56784x2X6] 

-167r^^[1020460x| + 177855x| + 1047914x2X4 + 19362x6] 

+1287r^^[220109x2^ + 70513x4] - 2558976077^^2 + 9265536 tt^^ . 

EM tensor & TC density: 

G'o\o) = :^^^, Gf{e) = 2 

GfiO) = 4[12x4 + x^] - 7r232x2 + 287r\ 

G^o\e) = A[2tItI - 6x|x4 + I5X3V4 - 40x2x2 - 25x2X3X5 - 625x2] 
+87r2[3x2^ - 11x2x1 + 73x|x4 + 6OX42 + 275x3X5] - 567r^[8x| + 5x1 + 40x2X4] 
+87r^[253x| + 27OX4] - 328077^2 + 1680 7r^° , 

Gf\e) = 4[AtItItI - 12x|x| + 3OX3VI - 80x2X4^ - I2X2VIX5 + 38x|x3X4X5 
-99x|x4X5 + 246X2X3X4V5 - I8X2V52 + I65X2X3V52 - 32OX2V4X52 - 45OX4VI 
+1125x3 x^ + 30x1x1x6 + 8IX3V6 - 8OX2V4X6 - 234x2x1x4X6 + 212x3X4X6 
-48x|x6 - 438x2x3X5X6 + 1530x3X4X5X6 - 9300x2X52x6 - 48X2V62 - 5022x32x62 
+10332x2X4X62 - 59940x6^] 
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-22t^T3T5 + 3T2Tln - 448r|r3r4r5 - 390r3r4V5 + 550t^t^ - 2325t^t^ 
+4225t2T4t| + 88x2^6 + 246T^r^rQ + AAAtIutq + I683T3V4T6 - 502OT2T4V6 
+II22OT2T3T5T6 + 22125x5^6 - 10854t|t6^ + 8730t4T6^] 

+87r^[42r2V| - 72r2r^ - lUr^n + 832r|T|T4 - 1804r|r| + 179x3^1 
-1696r2r| + 236r2V3r5 + 2943r|r5 - 3809r2r3r4r5 - 8435r2V52 
-13075r4r52 - 1952r2V6 - 10899r2r3V6 + 13288r2V4r6 + 5008r4V6 
-46035r3r5r6 + 38916r2r6^] 

+167r^[34T| - 480t|t| - 597t3^ + I79OT2V4 + 424T2T3V4 + 4665t|t| + 1182t| + 
328IT2T3T5 + 9435T3T4T5 + 17425t2t| + 2692t|t6 + I7979T3V6 
-2953OT2T4T6 - 25185x6^] 

+327r^[-540r2^ + 291x2^1 - 6338r2V4 - 2488r3V4 - 4222r2r4^ - 9908r2r3r5 
-9775T52 + 698r2V6 + 18322r4r6] 

+647r^°[2308T2^ + 1170t2t| + 9556t|t4 + 1012x4^ + 6345T3T5 - 1980t2T6] 
-167r^^[35992T| + 7683t| + 53924t2T4 + 372t6] 

+327r^^[36053r| + 15241t4] - 1148928 tt^ + 440128 tt^^ . 
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